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In curvaton-type models, observable non-gaussianity of the curvature perturbation would come 
from a contribution of the form {8a) 2 , where Sa is gaussian. I analyse this situation allowing 5a to 
be scale-dependent. The actual curvaton model is considered in more detail than before, including 
its cosmic uncertainty and anthropic status. The status of curvaton-type models after WMAP year 
three data is considered. 
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I. INTRODUCTION 

It is generally agreed that the primordial curvature 
perturbation £ is caused by the perturbation of one 
or more scalar fields, those perturbations being gener- 
ated on each scale at horizon exit during inflation. In 
curvaton-type models, a significant or dominant contri- 
bution to £ is generated after slow-roll inflation ends, by 
a field a whose potential is too flat to affect the inflation- 
ary dynamics. 

The body of this paper is in three sections. Section ITT1 
focuses on a perturbation of the form 



C(x) = Cinf(x)+Cr(x) 

C„(x) = bda (x) + Sa 2 (x) 



(1) 



(The last term is written with the compact notation 
Sa 2 = (Sa) 2 , which will be used consistently.) Follow- 
ing t ne spectrum, bispectrum and trispectrum of the 
perturbation are calculated, allowing for the first time 
spectral tilt in the spectrum of Sa. 

In Section IIIII the generation of ( is described using 
the SN formalism. In curvaton-type models £ is given by 
Eq. or its multi-field generalization. Non-gaussianity 
and scale-dependence are treated together, building on 
the separate discussions of 0, • 

The prediction for £ has cosmic uncertainty because it 
depends on the average value of the curvaton-like field 
in our part of the universe. One may assume that the 
probability distribution for this average typically is quite 
flat up to cutoff. The resulting probability distribution 
for £ (the 'prior' for anthropic considerations) is model- 
dependent. 

Section IIVI considers the actual curvaton model. A 
master formula for £ is presented, including all known 
versions of the model. Assuming that the curvaton con- 
tribution dominates and that the curvaton has negligi- 
ble evolution after inflation, the cosmic uncertainty and 
anthropic status of the curvaton model is described, ex- 
tending the recent work of Garriga and Vilenkin Q and 
Linde and Mukhanov @. Finally, in Section E| we look 
at the status of curvaton-type, in the light of the recent 
measurement of negative spectral tilt for the curvature 
perturbation. 

Standard material covered in for instance 0,0 is taken 
for granted throughout, with fuller explanation given for 



more recent developments. 

II. CALCULATING THE CORRELATORS 

For convenience it is assumed that the spatial average 
of Sa vanishes. 



Sa = 0. 



(2) 



This requirement is not essential, 1 because if Eq. were 
valid with some Sa =/= 0, one could arrive at Sa = by 
making the redefinitions 



Sa — > Sa — Sa 
b -> b + 2Sa. 

A. Working in a box 



(3) 



A generic cosmological perturbation, evaluated at 
some instant, will be denoted by by g(x) and its Fourier 
components by 



ffk 



x e 



: .k-x 



5(x). 



(4) 



The integral goes over a box of size L, within which the 
stochastic properties arc to be defined. To describe the 
cmb anisotropy the box should be much bigger than the 
size Hq 1 of the observable Universe. 

Since the box introduces periodic boundary conditions, 
one requires that physically significant wavelengths are 
much shorter than the box size, corresponding to kL S> I. 
One can then regard the wave vector k as a continuous 
variable. 

To describe the stochastic properties of cosmological 
perturbations within the box, one formally invokes an 
ensemble of universes and takes expectation values for 
observable quantities. The zero mode of each pertur- 
bation, corresponding to the spatial average within the 



In Q] and elsewhere, a contribution — Scr 2 was added to £ to make 
£ = 0. That has no effect on the calculations, which deal only 
with Fourier modes of f with nonzero wavenumber. 
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box, is not regarded as a stochastic variable. The nonzero 
modes have zero expectation value, (g^) = 0. (Both of 
these features are predicted by the inflationary cosmol- 
ogy.) It is usually supposed that the observable Universe 
corresponds to a typical member of the ensemble, so that 
the expectation values apply. 

Since the stochastic properties are supposed to be 
invariant under translations and rotations (reflecting, 
within the inflationary paradigm, the invariance of the 
vacuum) a sampling of the ensemble in a given region 
may be regarded as a sampling of different locations for 
that region. One can say then, that within the box of size 
L we are dealing with the actual Universe, and that the 
expectation values refer to the location of the observable 
Universe within the box. 

If \u(LHq) is not exponentially large, it should be safe 
to assume that our location within the box is typical. On 
the other hand, the observable Universe may be part of 
a very large region around us with the same stochastic 
properties; a region so large that Izi(LHq) can be ex- 
ponentially large. This is what happens within the in- 
flationary cosmology, if inflation lasted for an exponen- 
tially large number of Hubble times before our Universe 
left the horizon. If such a super-large box is used, one 
should bear in mind the possibility that our location is 
untypical. 

The spectrum P g (k) is defined by 

(5k5k') = (2 7 r) 3 ^(k + k')P 9 (fc). (5) 

It is useful to define V g = (fc 3 /2ir 2 )P g , also called the 
spectrum. After smoothing on a scale R, the variance is 

( 5 2 (x)) = jf* fr g (k). (6) 

The spectral index n g and the spectral tilt t g are defined 
as 



For constant tilt, V g oc /c' 9 . 

If V g (k) is sufficiently flat and the range of k is not too 
big, 

(9 2 )-P 9 - (8) 

If instead it rises steeply, (g 2 ) ~ T' g {R^ 1 ). In either case, 
the spectrum of a quantity is roughly its mean-square. 
This interpretation of the spectrum is implied in many 
discussions, including some in the present paper, but it 
should be applied with caution. 

On cosmological scales Vq is almost scale-invariant 
with -P c 1/2 = 5 x 1CT 5 . At the 2a level, the tilt is 
constrained by observation Q to something like = 
-0.03 ±0.04. 

If the two-point correlator is the only (connected) one, 
the probability distribution of <?(x) is gaussian. Non- 
Gaussianity is signaled by additional connected correla- 
tors. Data are at present consistent with the hypothesis 



that C is perfectly gaussian, but they might not be in the 
future. 

The bispectrum B g is defined by 

(3k l5 k 2 3k 3 ) = (2^) 3 ( 5( 3 )(k 1 +k 2 + k3) J B 9 (fc 1 ,fc 2 ,fc 3 ). (9) 

Instead of Bq it is more convenient to consider 
fahifa, k 2 ,k 3 ), defined by @ 2 

B ( {ki, k 2 , k 3 ) = ^/ NL [P c (fci)P c (*2) + cyclic] , (10) 
5 

where the permutations are of {hi, fc 2 , k 3 }. 
Current observation [l2l Il3l| gives at 2a level 

-27 < /nl < 121. (11) 

Absent a detection, observation will eventually [l4| bring 
this down to |/nl| ^ 1- At that level, the comparison of 
theory with observation will require second-order cosmo- 
logical perturbation theory, whose development is just 
beginning 

The trispectrum T g is defined in terms of the connected 
four-point correlator by as 

(3 kl 3k 2 3k35k 4 ) c = (2^) 3 <5 (3) (kj + k 2 + k 3 + k 4 )T 9 . (12) 

It is a function of six scalars, defining the quadrilateral 
formed by {ki, k 2 , k3, k4}. It is convenient to consider 
tnl defined by [lj 

T ( = ir NL P c (fc 1 )P c (fc 2 )P c (fc 14 ) + 23 perms. , (13) 

In this expression, ky = k^ + kj , and the permutations 
are of {ki, k 2 , k 3 , k 4 } giving actually 12 distinct terms. 

Current observation gives something like ^(| tnl ^ 
10 4 , and absent a detection PLANCK data will give 
something like tnl ^ 300. 



B. The correlators 

If £ is given by Eq. Q its spectrum, bispectrum and 
trispectrum are given by the following expressions, with 



2 The sign and the prcfactor make this definition coincide with the 
original one ll J in first-order cosmological perturbation theory, 
where /nl was defined with respect to the Bardeen potential 
which was taken to be $ = (In many theoretical works, in- 
cluding previous works by the present author, /nl is defined with 
the opposite sign.) At second order, which as we see later may be 
needed if |/nl| ^ 1, ^ and £ are completely different functions 
and /nl defined with respect to "3? has nothing to do with the 
/nl of the present paper. Unfortunately, both definitions are in 
use at the second-order level. 
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all higher connected correlators vanishing; 



p<* 
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^quad 
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23 perms. 

Tsa 2 
47T E 

IL- 1 

23 perms. 



(14) 
(15) 
(16) 

(17) 

(18) 

(19) 

(20) 
(21) 

(22) 



,3 Po{p)P„{pi)Pa{p2)Po{p2i) 



d A p 



p z p\p\p\^ 



(23) 



The 24 terms in Eq. (|22|l are actually 12 pairs of identical 
terms, and the 24 terms in Eq. (|23|l are actually 3 octu- 
plets of identical terms. In the integrals p\ = |p — ki|, 
pi = |p+k 2 | andp24 = |p+k 2 4|. The subscript indi- 
cates that the integrand is set equal to zero in a sphere of 
radius L _1 around each singularity. The integral (jl7|l was 
given in and the integrals H20|) and Q23JI were given 
in except that P a was taken to be scale- independent. 
The term ^ inear was given in [n) and the term T 1 - 1 j lear 
was given in 0] . 

In the language of field theory, these expressions are 
obtained by contracting pairs of fields, after using the 
convolution 



(2 



^ J d 3 qSa (s 5a k - cl . 



(24) 



The terms labeled 'quad' are generated purely by prod- 
ucts of three Sa 2 terms, while the terms labeled 'linear' 
are generated by a product of two bSa terms and one Sa 2 
term. Their evaluation is best done using Feynman-like 
graphs (cf. l2Qjl. The terms labeled 'linear' come 

from tree-level diagrams, while those labeled 'quad' come 
from closely-related one-loop diagrams. 

Since Sa is gaussian, its stochastic properties are de- 
termined entirely by its spectrum P a . The correlators 
(hence all stochastic properties) of da 2 are also deter- 
mined by P„. By examining the large-p behaviour of 
Eqs. I|17fl and (|20f> . we see that the correlators of Sa 2 
on a given scale are insensitive to the spectrum of Sa on 



much smaller scales. 3 In contrast with the case of quan- 
tum field theory, there is no divergence in the ultra-violet 
(large fc) regime. 

By examining the behaviour of Eqs. (|17l) and H2(J|) near 
the singularities, we see that with sufficiently large posi- 
tive tilt, the correlators of Sa 2 on a given scale are insen- 
sitive to the spectrum of Sa also on much bigger scales. 
More will be said later about this infra-red regime. 

The integral l|17(l can be evaluated exactly [TtJ to give 



P 5a 2{k) = 4P 2 ln(kL). 



(25) 



The integral H2U|) can be estimated as follows. Focus- 
ing on the singularity p = 0, one can consider a sphere 
around it with radius fc a bit less than min{fci, fc 2 }. The 
contribution from this sphere gives 



B. 



Sa 2 



327T 4 

prp 



dp 

L-i P 



32tt 4 

Ki^2 



ln(fcL). 



(26) 



A similar sphere around each of the other two singular- 
ities gives a similar contribution, and the contribution 
from these three spheres should be dominant because 
the integrand at large p goes like p~ 9 . (Applying this 
argument to the integral (|17|) happens to give exactly 
Eq. Evaluating /nl we arrive at the estimate pj 



3, 
5 

3 /.quad 



V 2 3 

u2 r a , u /•quad 

° Pi + 5 /nl 
P 3 

4^§ln(fcL) 

/ T~ 



P S g 

21n(fcL) V P( 



,-1/2 



(27) 
(28) 

(29) 



with k = min{/c 1 j_fc 2 ; ^3}- 
trispectrum gives [jj 



A similar estimate for the 



tnl 

quad 



.Pi 



quad 



16-fm(fci) 
' c 

1 (Via* 



lnfci V Pc, 



1 



(30) 
(31) 
(32) 



with k — minjfci, kj m }. 

It is easy to repeat these estimates for the case of con- 
stant nonzero tilt, t a = dh\P a / 'dlnfc. To avoid rather 
cumbersome expressions, I give the result for the bispec- 
trum only in the regime where the ki have an approxi- 
mate common value fc, and the result for the trispectrum 



3 To be precise, Eqs. 1171 . 121)1 . and 1231 converge at p 3> k if the 
tilt try is below respectively 3/2, 2 and 9/4. These conditions are 
well satisfied within the inflationary paradigm, which makes \t a \ 
well below 1. 
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only in the regime where both fcj and fcy have an approx- 
imate common value k. Then, the only effect of tilt is to 
replace ln(fcL) by 



y(kL) 



k dp /p\t° _ 1 - {kL)- l ° 
L -i p \k) 



(33) 



The following limits apply 



i (t a » l/ln(fcL)) 

y(kL) = { H^L) (\t a \ < l/ln(fcL)) . (34) 
(t a «-l/ln(fci) 

The expressions for the correlators in terms of y are; 
V Ca = b 2 V a +V Sa 2 (35) 



V 5a 2 ~ 4y(kL)Vt 

V Ca ~ (6 2 + 4yP CT ) 

3 . P» 
5' 



C 



, 1 / 1 



7> 2 2 V y(fcL) V y py 2 



c 

r NL ~ (6 2 + 4yV„) 

. 4^ ■ 1 ^ 



The tilt of C, is given by 



(36) 
(37) 

(38) 



(39) 
(40) 
(41) 

(42) 



with 



2t a (U > l/m(fc£)) 

= < l/]n(fci) (M « l/ln(fcL)) (43) 
t ff (t a « -l/ln(fcL) 

For zero or negative tilt, increasing the box size has an 
ever-increasing effect on the correlators. For positive tilt 
we can use a maximal box such that ln(A:L max ) l/t a 
and y = \jt a are good approximations, giving 



v 6a * c -vt 

5/NL - ^+2^ (^J 



7?3 

t nl = 4fe 2 — I + t a 



1 



(44) 
(45) 

(46) 



C. Working in a minimal box 

To minimize the cosmic uncertainty of the correlators, 
one might wish to choose the box size to be as small as 



possible, consistent with the condition LHq 3> 1 which 
is required so that it can describe the whole observable 
Universe 0|- How big LHo has to be depends on the 
accuracy required for the calculation of observables, us- 
ing the curvature perturbation as the initial condition. 
As the equations required for that calculation involve k 2 
rather than k it may be reasonable to suppose that very 
roughly 1% accuracy will be obtained with LHq ~ 10 and 
0.01% accuracy with LHo ~ 100. Even with the latter, 
the minimal box size L m i n corresponds only to 



ln(L min if ) - 5. 



(47) 



The range of cosmological scales is usually taken to be 
only Alnfc ~ 14, going from the size Hq 1 ~ 10 4 Mpc 
of the observable Universe, to the scale 10 -2 Mpc which 
encloses a mass of order 10 6 M Q and which corresponds 
to the first baryonic objects. 

Cosmological scales therefore correspond to roughly 



ln(fcL min ) - 5 to 20. 



(48) 



Let us see how things work out with the minimum box 
size. We saw that the dependence on the box size is 
through the function y(kL) given by Eqs. ll/i-'il) and 13411 . 
With the minimal box, this function is of order 1 on 
all cosmological scales, provided that < 1/20. This 
bound on t a is more or less demanded by observation if 
Cinf is negligible, but it can be far exceeded if £inf dom- 
inates. In the latter case I will allow only positive tilt, 
since strong negative tilt looks unlikely as seen in Section 
IIII Gl Then y is at most l/t a , and hence still roughly of 
order 1. 

To go further with the minimal box, it will be enough 
to consider the two extreme cases, that the correlators 
arc dominated by cither their 'linear' or their 'quad' con- 
tributions. From Eqs. (JSTJl, lEEJl, and the former 
case corresponds to 



Va^b 2 



(49) 



Given Eq. 10, this corresponds to the linear term of <^ a 
dominating, while the opposite case corresponds to the 
quadratic term of £ CT dominating. 

If the linear term dominates and £i n f is negligible, 



T c = b 2 V a 

-/nl = b- 2 
1 

r NL = (36/25)/ 2 L . 



(50) 
(51) 
(52) 



This is the usually-considered case. With a change of 
normalization one can write pT| (, = 5a + |/nl^ct 2 - 

If instead the quadratic term dominates it cannot 
dominate £ or there would be too much non-gaussianity. 
Indeed, given the interpretation JHJ), the non-gaussian 
fraction is 



Vs 



1/2 



1% 



1/2 



(53) 
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and 

r-n g ^(|/NL|^ 1/2 )^(TK L n)'- (54) 

The present bound |/nl| < 121 requires r ng < 0.2, but 
the present bound on tnl < 10 4 requires r ng < 0.07. 
We see that if the quadratic term dominates, the present 
bound on the trispectrum is a stronger constraint than 
the one on the bispectrum. Absent a detection, the post- 
COBE bounds on the bispectrum and trispectrum will 
lead to about the same constraint, r ng < 0.04. 

With the quadratic term dominating, the present 

1 /2 r 

bound on tnl gives < 3 x 10 . This, though, 

is on the fairly large cosmological scales probed by the 
cmb anisotropy. With positive tilt (A; max ) could be 
much bigger, even not far below 1 leading to black hole 
formation. 



are the spectrum, bispectrum and trispectrum defined 
within a particular small box of size L, we should have 

(P& 3 ) - V Sa , (57) 
(Bt 2 ) = B 5a2 (58) 
(T£ 3 ) = T s& >, (59) 

where the right hand sides and the expectation values 
refer to the big box. One can verify this explicitly using 
Eqs. 1)15(1 . 1)18(1 . and ((21(1 . Indeed, these equations apply 
to a box of any size. The hierarchy k ^> L _1 3> M _1 al- 
lows one to evaluate the changes in the integrals induced 
by the change M — * L, and to verify (for any form of the 
spectrum) that this change is precisely compensated by 
the change b 2 — > {b 2 L ). 

III. THE INFLATIONARY PREDICTION 



D. Running the box size 

We have found that the stochastic properties depend 
on the size of the box in which Eqs. and J2J are sup- 
posed to hold. This seems to be incompatible with a 
basic tenet of physics concerning the use of Fourier se- 
ries, that the box size should be irrelevant if it is much 
bigger than the scale of interest. 

This situation was discussed in on the assumption 
that V a is scale-independent, and it is easy to extend the 
discussion to the case of an arbitrary V a (k). The crucial 
point is that 8a is supposed to vanish within the cho- 
sen box of size L. Let us imagine now that this box is 
within a much bigger box of fixed size M, and see how 
things vary if the size and location of the smaller box 
are allowed to vary. We have in mind that the small box 
will be a minimal one, and in the case of constant posi- 
tive tilt the big box might be the maximal one satisfying 
\\i(MHa)t a 3> 1. Defined in the big box, £ has the form 
(JTJ, with some coefficient b and with 8a = 0. 

Focus first on a particular box with size L, and denote 
quantities evaluated inside this box by a subscript L. In 
general 8a l will not vanish, and absorbing its expectation 
value into b using Eq. © we find 



b L =b + 25a L . 



(55) 



Now, instead of considering a particular small box, let its 
location vary so that 5a l becomes the original perturba- 
tion 8a smoothed on the scale L. Then 



(bl) 



= b 2 + A{8a 2 ! 



{8a 2 



dk 



M- 1 



Va{k), 



(56) 



where the expectation values refer to the big box. 

The operations of smoothing and taking the expecta- 
tion value commute. Therefore, if V 



5u 2 ' 



B% 2 and Tf ■ 

ba A da- 



In this section we see how Eq. (QJ and its generaliza- 
tions may be predicted by inflation. We take the rele- 
vant scalar fields to be canonically normalized, and focus 
mostly on slow-roll inflation with Einstein gravity. The 
latter restriction is not very severe because a wide class 
of non-Einstein gravity theories can be transformed to an 
'Einstein frame' 0. 



A. The curvature perturbation 

The curvature perturbation £ is defined on the space- 
time slicing of uniform energy density p through the met- 
ric [lElElliiliil 



= a 2 (t) 



,2£(x,t) 



7»j(x), 



(60) 



where 7 has unit determinant. Within the inflationary 
cosmology, if the tensor perturbation is negligible, 7^ — 

By virtue of the separate universe assumption, the 
threads of spacetime orthogonal to the uniform density 
slicing are practically comoving. As a result, a comoving 
volume element V is proportional to a 3 (x) where 



i(x,t) = a(i)e C(x '* } . 



(61) 



This means that a(x, t) is a locally-defined scale factor. 

An alternative definition of Q refers to the slicing where 
the metric has the form gij = a 2 ^ij (with again 7 having 
unit determinant). This is usually called the spatially flat 
slicing, which it is if the tensor perturbation is negligible. 
Linear cosmological perturbation theory gives £ in terms 
of Sp on the spatially flat slicing, 



c 



-H 



Sp 



3(P + P) 



(62) 



If /nl turns out to be of order 1 though, it will be neces- 
sary to go to second order and then the non-perturbative 
definition 1(60(1 becomes more useful 0,H3- 
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The linear expression <|62(l is convenient if the fluid is 
the sum of fluids, each with its own P(p). Defining on 
flat slices the constants 



one has 



C(*) = 



8 Pi 



3(p» + PO ' 



E(Pi+_P)6 
(p+p) 



(63) 



(64) 



By virtue of the separate universe assumption, the 
change in the energy pV within a given comoving vol- 
ume element is equal to —PdV with P is the pressure. 
This is equivalent to the continuity equation 



q a(x,t) 

P = -3 / Ts (P + P) 

a(x, t) 

d(x, t) a ' 
— - — - — — = ~ + (. 

a(x, t) a 



(65) 
(66) 



Remembering that this holds on uniform-densit y sl ices, 
we see that C is conserved [H Q El El HI M HI 
during any era when P is a unique function of p. That 
is guaranteed during any era where there is practically 
complete radiation domination (P = p/3), matter domi- 
nation (P = 0) or kination (P = p). The generation of £ 
may take place during any other era. 



B. Generating the curvature perturbation 

The idea is that during inflation, the vacuum fluctua- 
tion of each light field becomes, a few Hubble times after 
horizon exit, a classical perturbation. (To keep the lan- 
guage simple I shall loosely say that the classical pertur- 
bation is present at horizon exit.) The observed curvature 
perturbation, present a few Hubble times before cosmo- 
logical scales start to enter the horizon, is generated by 
one or more of these classical field perturbations. 

Opportunities for generating £ occur during any era 
when there is no relation P(p). Generation was origi- 
nally assumed to take place promptl y at horizon exit in 
a single-component inflation model |2lL I22I |23| ]. Then 
it was realised p4| that in a multi-component inflation 
model there will be continuous generation during infla- 
tion, but it was still assumed that the curvature pertur- 
bation achieves its final value by the end of inflation. 

The term 'curvaton-type models' in this paper denotes 
models in which a significant contribution to the curva- 
ture perturbation is generated after the end of slow-roll 
inflation, by the perturbation in a field which has a negli- 
gible effect on inflation. If the curvaton-like contribution 
is completely dominant then the mechanism of inflation 
is irrelevant. 

The curvaton model itself was the original proposal. In 
this model, the oscillating curvaton field leads to a second 
reheating, and the curvature perturbation is caused by 



erated by an inhomogeneity in any or all of the deca; 
rate 0, |35|, the mass j3|| or the interaction rate 



the perturbation in the curvaton field |25ll26ll27ll2al29||. 4 

Alternatives to the curvaton model, which still use a 
reheating, are to have the curvature perturbation gen- 
cay 

m 

of the particles responsible for the reheating. In that 
case the reheating can be the first one (caused by the 
scalar field(s) responsible for the energy density during 
inflation) or alternatively the particle species causing the 
reheating can be a fermion |38j . Other opportunities for 
generating the curvature perturbation occur at the end 
of inflation [32|. during preheating [33j, and at a phase 
transition producing cosmic strings [39 ]. 

To describe the generation of the curvature perturba- 
tion I make two assumptions. First, that the evolution 
of perturbations on a given cosmological scale k can be 
described using an idealized universe, which is smooth on 
some scale a bit smaller than 1/fc. Second, that the local 
evolution of this idealized universe in the super-horizon 
regime can be taken to be that of some unperturbed uni- 
verse. 

The first assumption is routinely made in cosmology, 
on both super- and sub-horizon scales. The second as- 
sumption is the separate universe assumption p3ll40ll4l| . 
which amounts to the statement that the smoothed uni- 
verse becomes locally isotropic and homogeneous when 
the smoothing scale is much bigger than the horizon. 
Given the content of the Universe at a particular epoch, 
it may be checked using cosmological perturbation the- 
ory, or else using the gradient expansion El El El D3 
with the additional assumption of local isotropy. Local 
isotropy is more or less [43 guaranteed E|| by inflation. 
Independently of particular considerations, the separate 
universe assumption has to be valid on a scale a bit larger 
than Hq 1 , or the concept of an unperturbed FLRW Uni- 
verse would make no sense. It will therefore be valid on 
all cosmological scales provided that all relevant scales in 
the early Universe are much smaller, which will usually 
be the case |47|. 

I assume slow-roll inflation, noting possible generaliza- 
tions along the way. The light fields fa are defined as 



4 Earlier papers [30] considered the generic scenario, in which a 
light scalar field gives a negligible contribution to the energy den- 
sity and the curvature perturbation during inflation, but a signif- 
icant one at an unspecified later epoch. Such a scenario becomes 
a curvaton-type model if that epoch is before cosmological scales 
start to leave the horizon, otherwise it may be an axion-type 
model giving a cdm isocurvature perturbation. Among curvaton- 
type models, the curvaton model is the one which generates the 
curvature perturbation from the perturbation in the am plitude 
of the oscillating curvaton field. It was described in l25l . a nd a 
for mu la equivalent to the estimate £ ~ 8a jo was given in 12(1 . 
In l27l the curvaton model was advocated as the dominant cause 
of the curvature perturbation and a precise c alcu lation of 
was made allowing for significant radiation. In 28] a significant 
inflaton cont ribu tion was allowed. The first calculation of /nl 
was given in l2?l . The curvaton mech anism with a pre- big- bang 
instead of inflation was worked out in 1 .3 1 1 . 
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those which satisfy flatness conditions; 

£i<l |%-|<1, 

where 



M 2 ^ 



with Vi = dV/dfa and V tJ = d 2 V/ 'dfcdfa . 
The exact field equation for each light field, 

4n + + Vi=Q, 

is supposed to be well-approximated by 

3Hfr = -Vi. 



(67) 

(68) 
(69) 

(70) 
(71) 



In these expressions H is the Hubble parameter, re- 
lated to V by 



3M|if 2 



V- 



(72) 



By virtue of Eqs. JBSJ, and lf7I|l. this becomes 

3MpH 2 ~ V, and i/ is slowly-varying corresponding to 
almost-exponential inflation; 

1 (iff 

2e < 1, (73) 



if 2 dt 



with 



(74) 



(Despite the notation, it is y 7 ^ and not which trans- 
forms as a vector in field space.) If inflation is almost 
exponential but not necessarily slow-roll, it may still be 
useful to define the light fields by Eqs. f^-tp, with V 
in the denominator replaced by 3MpH 2 . 

Focusing on a given epoch during inflation, it may be 
convenient to choose the field basis so that one field <f> 
points along the inflationary trajectory. I will call it the 
inflaton, which coincides with the standard terminology 
in the case of single-component inflation. I will denote 
the orthogonal light fields (assuming that they exist) by 
<7j. In an obvious notation e CTi = and e — initially. 
One may define also rj = r\^. 

From Eqs. an d CD; the gradient of the po- 

tential is slowly varying; 

The possible inflationary trajectories are the lines of 
steepest descent of the potential. The trajectories may be 
practically straight (single-component inflation) or sig- 
nificantly curved in the subspace of two or more light 
fields (multi-component inflation, called double inflation 
in the case of two fields). The terms single- and multi- 
component refer to the viewpoint that the inflaton field 
is a vector in field space. For single-component inflation, 
e ai «e, and the inflaton field <f> hardly changes direction 
in field space , so that one can choose a practically fixed 
basis {(/), <Ji}. 



C. The 5N formula 

To evaluate the curvature perturbation generated by 
the vacuum fluctuations of the light fields, we can use 
the SN formalism 0, H El It gives C(x,t), in 
terms of the light field perturbations defined on a flat 
slice at some fixed 'initial' epoch during inflation; 



(76) 



Keeping terms which are linear and quadratic in Sipi, the 
time-dependent curvature perturbation is 



C(x,i) = Mf(&(x),p(t)) 



(77) 



Here, N(<f>i,p) is the number of e- folds, evaluated in an 
unperturbed universe, from the initial epoch to an epoch 
when the energy density p has a specified value. In the 
second line, Ni = dN/dfa and = d 2 N/d(j)id(pj, both 
evaluated on the unperturbed trajectory. In known cases 
the first two terms of this expansion in the field pertur- 
bations are enough. 

The curvature perturbation £(x, t) is independent of 
the 'initial' epoch. At that epoch let us work in a basis 
(4>,(7i) so that Ni — (JV^, iVo-j). The contribution of the 
inflaton to £ is time-independent because S(f> just corre- 
sponds to a shift along the inflaton trajectory. It is given 
by 



Cinf 

N H /Nl 



1 V 
77 - 2e. 



(78) 
(79) 
(80) 



Redefining (f> we learn that the non-gaussianity coming 
from 84> 2 is negligible (|/nl| <C l)- 5 We can therefore 
write 



(81) 



The fields correspond to shifts orthogonal to the 
inflaton trajectory. The perturbations Sai give initially 
no contribution to the curvature perturbation; in other 
words the <Ji are initially isocurvature perturbations. If 
the N (Ji were evaluated at a final epoch just a few Hub- 
ble times after the initial one, they would be practically 
zero because the possible trajectories can be taken to be 
practically straight and parallel over such a short time. 
The crucial point is that the N ai can subsequently grow. 



The three-point correlator of 6<f> also I20U55I gives |/nl| C 1. A 
different proof that single-component inflation gives |/nl| <C 1 
was given earlier by Maldacena 0. 
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Such growth has nothing to do with the situation at hori- 
zon exit, in particular there is no expression for the N ai 
in terms of the potential gradient analogous to Eqs. 179(1 
and ijHUjl . 

In any case, £ settles down at some stage to a final 
time-independent value, which persists until cosmological 
scales start to enter the horizon and is constrained by 
observation. This value is the one that we studied in the 
previous section, focusing on Eq. which is obviously 
a special case of Eq. <j77|) . 

Although we focused on slow-roll inflation, the basic 
formula (|77|l can describe the generation of £ from the 
vacuum fluctuation for any model of almost-exponential 
inflation, which need have nothing to do with scalar fields 
or Einstein gravity. The only requirement is that the 
specified light fields fa, satisfying Eqs. ((67(l - l(tj9J) with V 
in the denominator replaced by 3MpH 2 , determine the 
local evolution of the energy density and pressure until 
the approach of horizon entry. Going further, the fa need 
not be the light fields themselves, but functions of them, 
evaluated at an 'initial' epoch which might be after the 
end of inflation. 



to give dEl 

/ r lnm N n N m 2 \ 

C V NiNi MpNiNi J t V ' 

with the right hand side evaluated at horizon exit. 

At horizon exit, let us use the basis fa Oi. If the infla- 
ton contribution 5(f) dominates one recovers the standard 
predictions; 

p <^4^© 2 (87) 

k - ** = C-n - 6e), ■ (88) 

In this case, the tensor fraction Vt/V^ is r = 16e. 

The contributions <J(7j to Vq are positive making [(| 
r < 16e. If one or more of them completely dominates, 
r will be too small to ever observe. If just one of them 
dominates (call it a) dominates, 

n- = K. (f ) 2 m 

t { ~ („ = (2,„ - 2e), . (90) 



D. The gaussian approximation 



E. Non-gaussianity in curvaton-type models 



In the limit where the potential of the light fields is per- 
fectly flat, inflation is exponential (de Sitter spacctimc) 
with constant H . At least up to second order in cosmo- 
logical perturbation theory |55| , the perturbations in the 
light fields arc then generated from the vacuum without 
back-reaction. Each perturbation is classical, and a few 
Hubble after horizon exit is practically time-independent 
with the spectrum [l0( (H/2ir) 2 . Keeping a slow vari- 
ation of H and finite values of the parameters and 
rjij , the three-point correlator of the light field perturba- 
tions has been calculated by Seery and Lidsey [5j|, and 
its effect has been shown indeed to be negligible in some 
particular cases [2(], H^, |5(| . In the following we take the 
light fields to be gaussian at horizon exit. 

Since (, is almost gaussian it will presumably be a good 
approximation to keep only the first term of Eq. (|77|l . Let 
us recall briefly the situation for that case. Focusing on 
a particular scale, we can take the initial epoch to be the 
epoch of horizon crossing. (The horizon-crossing trick.) 
Then 



C(k)=^iV«^(k), 



(82) 



where a star denotes the epoch of horizon exit. Using the 
summation convention. 



V c {k) = (H,/2tt) 2 N„N u 



(83) 



To evaluate the tilt from this equation one may use the 
slow-roll expressions 



d\nk 



Ml 
V 



-Vi- 



d 



NiVi = MpV, 



(84) 
(85) 



Now I consider the quadratic term of Eq. 1(77(1 , allow- 
ing for the first time the possibility of significant scale 
dependence. 

The horizon-crossing trick for evaluating scale depen- 
dence may not work when the quadratic term is in- 
cluded, because the convolution 1(24(1 involves a range of 
wavenumbers. It can still be used if the correlators of 
interest involve only a narrow range of scales ki ~ fc, and 
one just wishes to evaluate the dependence on the overall 
scale k. Otherwise, one should take the 'initial' epoch for 
Eq. ((77(1 to be a fixed one, after all relevant scales have 
left the horizon. 

To handle this situation I suppose that inflation is 
single-component, and that there is just one non-inflaton 
field a which is going to be the curvaton-like field. (With 
small modification the discussion still applies to multi- 
component inflation models generating negligible non- 
gaussianity.) Then 



C = Cinf + Ca 
Cinf — N^Stj) 



(a 



NJa + -N aa 5<j 



(91) 
(92) 

(93) 



Up to the normalization of <5<r, Eq. I(93|) is the same as 
Eq. <m>. 

Since e CT <^ e, the perturbations on flat slices satisfy 

M 



i 

H dt 
1 dSa 
H~df 



- (rj - 2e) 5cf) 

-VcrcrSa, 



(94) 
(95) 
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The second equation is just the first-order perturbation 
of the unperturbed slow-roll equation 3Ha = — V a . That 
equation therefore applies locally. 
The solutions are 



dt i 



8(7 



,dl 



8(7* 



(96) 
(97) 



These are to be evaluated at a fixed t which will be the 
'initial' epoch for use in Eqs. (|92|l and l|93l) . Keeping only 
the linear terms and evaluating the spectrum, we recover 
Eqs. (|88J) and l|9()(l for the spectral indices. 

From now on the focus will be on the quadratic poten- 
tial 



V(a) 



1 



Also H is supposed to be sufficiently slowly varying that 



3#r 



(99) 



where H* is now the practically constant value of H dur- 
ing inflation without reference to a particular epoch. 

As with any situation involving scalar fields in the early 
Universe, one has to remember that the effective poten- 
tial V(a) can be affected by the values of other scalar 
fields and change with time. In particular, supergrav- 
ity gives for a generic field during inflation \r\ aa \ ~ 1. 
This marginally violates the slow-roll condition and cor- 
responds roughly to to* ~ iJ* and t„ ~ 1. The field 
which dominates ( must have \n aa \ < 0.01. In any case, 
the mass to* appearing in Eq. (|98|l will generally not be 
the true mass to, defined in the vacuum. 



F. Maximum wavenumber and the 'initial' epoch 

The classical curvature perturbation £ is generated up 
to some maximum wavenumber fc max . This maximum is 
generally taken to correspond to a scale 1 / fc max far below 
the shortest scale of direct cosmological interest discussed 
in Section Til CI Nevertheless the value of fc max matters. 
It represents the shortest possible scale for the formation 
of primordial black holes 6 and the shortest scale on which 
matter density perturbations can exist. As we shall see, 
short-scale perturbations in the curvaton density con- 
tribute to its mean density, and hence indirectly to the 
magnitude of £ which the curvaton model generates. 

If £ is created during inflation, /c max is the scale k e 
leaving the horizon at the end of inflation. It is given by 



k c 
Ho 



JV 



(100) 



6 See for instance references in 
k > fc max is also considered. 



where the quantum regime 



where N is the number of e-folds of slow-roll inflation af- 
ter the observable Universe with present size Hq 1 leaves 
the horizon. For a high inflation scale with continuous 
radiation domination afterward N ~ 60. To make fc" 1 
comparable with the shortest cosmological would require 
N ~ 14 (see the discussion after Eq. Q47[) ) which is hard 
to achieve. For this reason it seems to have been as- 
sumed in all previous discussions that fc max will be far 
below cosmological scales. 

That assumption is not justified if C is created by a 
curvaton-type mechanism long after inflation. After in- 
flation, the perturbation 8a in the curvaton-type field 
redshifts away on scales entering the horizon. Therefore, 
fcmax is the scale entering the horizon when Q is created. 
(See |l7^ for the same phenomenon in the axion case) 
This scale leaves the horizon long before the end of in- 
flation. To handle that situation, I will still equate fc max 
with k e given by Eq. I|100|) . but define N as the number 
of e-folds of relevant inflation after the observable Uni- 
verse leaves the horizon, 'relevant' meaning e-folds which 
produce perturbations on scales bigger than fc max • De- 
manding only that fc max is below the shortest cosmolog- 
ical scale we can allow a range 



14 < N < 60. 



(101) 



Taking the extreme values N = 60 and (say) t a = 0.4 
gives e m " ~ 10 11 . 

The end of relevant inflation is the appropriate 'initial' 
epoch for use in Eq. (|93|l . The spectrum of a e is 



k 



2tt ) \kj \2tt 

(102) 

The factor e Nt " is of order 1 if Nt G < 1. This is more 
or less demanded by the observational bound on t$ if 
the curvature perturbation is dominated by the curvaton 
contribution. In the opposite case large tilt is allowed, 
making e Nta exponentially large as we saw earlier. 
The unperturbed field at the end of inflation is 



N L = ln{k e L), 



(103) 
(104) 



where a is the practically unperturbed value of the cur- 
vaton field within the box of size L when it leaves the 
horizon, and is the number of e-folds of relevant in- 
flation after the box leaves the horizon. 



G. Cosmic uncertainty 

If inflation lasts for enough e-folds before the observ- 
able Universe leaves the horizon, the stochastic formalism 
|59j allows one to calculate the probability distribution of 
a, for a random location of our Universe. If our location 
is typical the actual value of a 2 will be roughly (a 2 ). 
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If the variation of H is so slow that it can be ignored, 
one arrives at a particularly simple probability distribu- 
tion which is described in this subsection. Generaliza- 
tions to allow for varying H are given for instance in 



1. The Bunch-Davies case 

In the case of constant positive tilt the result can be 
obtained from the formalism already presented. To do 
this, one works in a maximal box and assumes that the 
average value of a rolls down to a practically zero value 
well before the observable Universe leaves the horizon. At 
any subsequent epoch, the local value er(x) has a gaussian 
probability distribution with variance 



(^(x)> 



2^ 

H,, 
2^ 



2 dk f k ' 

k \ aH* 



1 



87t 2 to 2 



(105) 
(106) 



This is the case considered by Bunch and Davies 10, 21] . 

Working within a smaller box with size L (thought 
of as a minimal one), c(x) has a spatial average and a 
perturbation; 



<t(x) = a L + <5ctl(x), 



(107) 



where the classical perturbation includes all wavenum- 
bers k < aH. The mean square within that box is 



(Sat 
aH * dk r 



(108) 
(109) 



For a random location of the small box (within the max- 
imal box) each term of Eq. I|107|) has a gaussian distri- 
bution. Adding the two variances gives 



iH, 



dk 
~k 



(110) 



which agrees with Eq. (|1()5|> . 

When the minimal box first leaves the horizon the per- 
turbation 8a is negligible. For a random location of the 
minimal box, the variance of the unperturbed value a 
is then practically equal to the Bunch-Davies expression 

We were defining the curvature perturbation i|93|) 
within a minimal box, because that has general appli- 
cability. In the Bunch-Davies case we can instead use a 
maximal box, big enough to ensure the condition a = 
before the observable Universe leaves the horizon. In that 
case, N„ may vanish leaving only the quadratic term of 
Eq. (|93[1 . This will happen if a — ► — a is a symmetry of 
the theory, and it happens anyway in the actual curvaton 



model because N a is then determined directly by the po- 
tential. Then the correlators are given by Eqs. (|44|) (|46[l . 
with 6 = and V$o-2 — V^ a ; 



- 1+1/2 



,1/2 



tnl — to 



(111) 
(112) 

(113) 



2. The general case 



In general, Eq. H98|) will contain higher terms 



d>A IU P 



(114) 



They will not affect the Bunch-Davies result provided 
that 



a « n 

d-4 

« tl 



(115) 
(116) 



For an arbitrary potential V(o~), assuming still that in- 
flation with practically constant H lasts for long enough, 
the probability of finding a in a given interval is |59l l60| 



P(a)da oc exp(-87r 2 U/3iJ^)dcr. 



(117) 



The potential is V(a) is evaluated with all other rele- 
vant fields fixed, with the convention that its minimum 
vanishes. 

There are two simple cases. If V(o~) increases until 
it becomes at least of order H%, the probability distri- 
bution is more or less flat out to a value cr max , such 
that V ~ 3H^/8ir 2 . 7 For the Bunch-Davies case one can 
take cr^a^ to be the variance (i?*/27r) 2 /t cr of the gaussian 
probability distribution. In this case the tilt is positive 
and it can be strong, corresponding to the lightness con- 
dition t a <C 1 being only marginally satisfied. 

Instead, a might be a PNGB with the potential 



V(a) = A 4 cos 2 (7r(j//). 



(118) 



If A <C the probability distribution for a is extremely 
flat within the fundamental interval < a < f. The 
effective mass at the maximum and minimum of the po- 
tential is rn* = 7rA 2 // giving the maximum tilt as 



t„ 



A!. A 2 
i^T 2 ' 



(119) 



7 This d iffers from the estimate of the typical value given in Ifill 

MM 



11 



To have a reasonable probability for being at the maxi- 
mum requires A -C -ff* , and to have weak self-coupling so 
that the semi-classical theory used here makes sense re- 
quires A -C /• Judging by this example, strong negative 
tilt looks unlikely and is not considered in the present 
paper. 

On the other hand, slight negative tilt consistent with 
observation is possible. The probabilities for being at the 
maximum and minimum of the potential are related by 



Pn 



exp 



8tt 2 A 4 



(120) 



and one can have say t a = —0.05 while keeping this ratio 
not too far below 1. 

It must be emphasized that the probability distribu- 
tion Eq. (|117fl is attained only if the variation of H is 
negligible, on the timescale for the rolling down of a to- 
wards its zero value. Whatever it is, the late-time prob- 
ability distribution is not relevant for the inflaton in a 
single-component inflation model, since its value a given 
number of e-folds before the end of inflation is obtained 
by integrating the trajectory 3H(j) = — V' . 

If eternal inflation takes place around some maximum 
of the potential, H will be practically constant during the 
eternal inflation and all light fields will attain the proba- 
bility distribution (|117fl except for the one driving eternal 
inflation. When eternal inflation ends and slow-roll in- 
flation begins, the fields orthogonal to the inflationary 
trajectory will have this probability distribution. In a 
single-component inflation model it will still apply when 
the observable Universe leaves the horizon, if the value 
of V then is not much lower than it was during eternal 
inflation. That may be the case if the potential has a 
suitable maximum, which is more likely than one might 
think The distribution l|117|) might also be attained 
P| if eternal inflation occurs high up on the chaotic in- 
flation potential V cx (f> 2 . A further possibility, not yet 
investigated, is that the distribution l|117|) is attained if 
eternal inflation occurs at a maximum of the potential 
with high V; possibly well-motivated I66ll67l realizations 
of that case would be Natural Inflation |68W along with its 
hybrid |69ll70j | and multi-component |7l| generalizations. 



on V(a) will have become negligible by this time. That 
will be more or less true if the effective mass up to that 
time has been < H(t). 

When the harmonic oscillation begins, p a is supposed 
to be negligible compared with the total p = p a + pi a {. 
The component p ln f (defined as the difference between 
p and pa) roughly speaking originates from the decay of 
the inflaton but there is no need of that interpretation. 
The curvaton contribution to the curvature perturbation 
is at this stage supposed be negligible. 

Eventually the harmonic oscillation decays. During at 
least some of the oscillation era, p m i is supposed to be 
radiation-dominated so that Pajpmi grows and with it 

Originally p; n f was supposed to be radiation- 
dominated during the whole oscillation era, but the 
model is not essentially altered if p ln f contains a sig- 
nificant contribution from matter. This matter might 
be the homogeneously oscillating inflaton field which 
decays only after the onset of the curvaton oscillation 
[EH E3, Elf, non-relativistic curvaton particles other 
non-relativistic particles which decay before the curvaton 
or any combination of these. 

The lightness of the curvaton field can be ensured by 
taking it to be a PNGB with the potential fTTBjl . This 
mechanism can work whether or not there is supersym- 
metry, and is easier to implement for the curvaton than 
for the inflaton [H EH, [zl . 

Several curvaton candidates exist which were proposed 
already for other reasons. Using such a candidate, one 
might connect the origin of the curvature perturbation 
with particle physics beyond the Standard Model, or even 
with observations at colliders and detectors. Among the 
candidates are a rig ht-handed sneutrino [lH Ell 151 

ta modulus |2£t l79j| (which might be a string axion 
lU), a Peccei-Quinn field [sUsl and an MSSM flat 
direction [79ll8^ |. The right-handed sneutrino possibility 
was actually discovered serendipitously [zl| by authors 
who were unaware of the curvaton model, which shows 
that the curvaton model is not particularly contrived. 



B. The master formula 



IV. THE CURVATON MODEL 
A. The setup 

The curvaton model [H |H HI |28| ^ alS ° 0) is a 
particular realization of Eq. I|93|l . The curvaton field a at 
some stage is oscillating harmonically about a = under 
the influence of a quadratic potential V = ^m 2 (j) 2 , with 
energy density p a cx a~ 3 . This stage begins at roughly 
the epoch 

H~m. (121) 

This mass m in these equations is taken to be the true 
vacuum mass, the idea being that the effect of other fields 



In this subsection the basic approach is that of |27l 
[29| . which works with the the first-order perturbation 
theory expression 1)62(1 . This is applied after the onset 
of the harmonic oscillation of the curvaton, when the 
cosmic fluid has two components. The final value of ( 
is taken to be the one evaluated just before the curvaton 
decays, which is taken to occur instantaneously on a slice 



Equations are derived on the assumption that these quantities 
are initially negligible. Presumably those same equations will 
provide a crude approximation even if the gr owth is negligible, 
due either to the curvaton decaying promptly l72t 17.*^ or to pi n f 
containing a negligible radiation component. 
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of uniform energy density, at an epoch H ~ T where F is 
the decay width. 



Keeping this basic setup, the treatment of |27],|29j will 
be generalized to allow for several possible effects. 9 The 
inflaton component p m f is not required to be purely ra- 
diation. The contribution of the inflaton perturbation to 
£ is not required to be negligible. The tilt t a is taken 
into account. Attention will focus on constant tilt which 
is either small (|i CT | < 10 -2 , or else large and positive, 
which is allowed if the curvaton contribution to £ is sub- 
dominant. In that case strong non-gaussianity will also 
be allowed. 

Evolution of the curvaton field after the end of inflation 
will be allowed. One possibility [64| for such evolution is 
the large effective mass-squared V aa ~ ±H 2 (t) predicted 
by supergravity for a generic field during matter dom- 
ination (though not |90j | during radiation domination). 
A more drastic possibility is for a to be a PNGB corre- 
sponding to the angular part of a complex field, whose 
radial part varies strongly [7^, 0, |^, [93 ■ 

In the presence of evolution, the oscillation may ini- 
tially be anharmonic, but after a few Hubble times the 
amplitude presumably will have decreased sufficiently 
that the oscillation is harmonic, making Eq. I|121|) an 
adequate approximation. (A detailed discussion for the 
analogous axion case is given in [13 •) 

The evolution is given by 



a + 3H(t)a + V„ = 0. 



(122) 



Since the inflaton perturbation just corresponds to a shift 
in time, the 'separate universes' are practically identical 
until after the onset of the oscillation, and Eq. (|122|l holds 
locally at each position. Let us define the amplitude 
CTos(x) at the start of the harmonic oscillation on a space- 
time slice of uniform energy density. Then a os (a e (x)) is 
a function only of a c , the x-dependence coming purely 
from the fact that a c (x) is not defined on such a slice. If 
V is quadratic (with a constant or slowly-varying mass) 
Cos(ce) is practically linear. 



9 Each effect was considered before, usually without any of th e 
others. Strong tilt and non-gaussianity were considered in Fl l2(l . 
The possible contribution of curvaton particles to p; n f was taken 
into account qu alitativ ely in |3- The Bunch-Davies case was 
considered in |26l |78jl. Cu rvaton evolution was partially taken 
into account in l2a l6l 153 (see also for a treatment using 
the SN formalis m)) . The po ssible contribution of £; n f was taken 
into account in 28, 61, 62]. All of this is at first order. The 
calculation t o secon d order in cosmological perturbation theory 
was done in l87t l88l (re-derived in using the SN formalism.) 
Also, the sudden-decay approximation was removed in |89j . at 
first order only. 



Knowing cr s(x) we can calculate 10 



P<x(x) 



— (c os + <5cr os (x)) 2 



Pa 



a, 



{8°Z 



Sp a = 



—— (2a 0S Sa 



Sa 2 



(123) 

(124) 
(125) 
(126) 



To go further we expand CT os (<7 e ) to second order in Sa e 
giving [H 



(5er os (cr c (x)) = a' oa 6a e (x) + -trl^(x), 



and 



Spa 

Pa 

q 
u 

V 



9 2 



mrp 2 



2q— + u 



qS<r c 



fcOog/cTos 

1 + <T os a os /a os 



(127) 



(128) 

(129) 

(130) 
(131) 

(132) 



If da has negligible evolution, or if a os (a c ) is linear cor- 
responding to a quadratic potential, then q = u = 1. It 
has been shown 93] that even slight anharmonicity could 
in certain cases give |u| -C 1, and as we have seen strong 
evolution is also possible making both q and u very dif- 
ferent from 1. 

Except in Section IIV Fl all calculations will be done 
with the minimal box size, corresponding to ln(L_ffo) not 
too far above 1. We need {5a 2 



k ° dk 



2^ 



y{Lk c 



This gives 



p = 1 + uq' 



\ 2ira J 



y(Lk e )e 



N L t a 



with 



y ~ mm 



N 

1/to 



(133) 
(134) 

(135) 
(136) 



A crude but usually adequate approximation is (<5<r 2 ) ~ 
H 2 . (See 01 f° r a similar estimate of the axion pertur- 
bation.) 



Recall (footnote 2) that there is no need to demand Sp a = 0. 
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We have been evaluating p a and its perturbation at the 
beginning of the oscillation, on a slice where p — p ul { is 
uniform. The curvature perturbation is given by Eq. 1641) 
in terms of the perturbations of the two fluids evaluated 
on the flat slicing. The curvaton perturbation on the 
uniform-density slicing is 



5 Pa _ ( 5p a 5p inl 

Pa 



Pinf + -Pinf 



(137) 



flat 



Each term in this expression is time-independent |29| . 
(Remember that the two-fluid description is only valid 
during the harmonic oscillation.) Using it, Eq. I)64|) eval- 
uated just before the curvaton decay gives 11 



C = Cinf+C 
Ca = f b - P ° 



3/ = 



Po 
Pa 



p + P p a + Pinf + Pinf P 



(138) 
(139) 

= n a . (mo) 



The approximation is adequate, because the 'sudden- 
decay' approximation has generally a significant error [89l | 
in the regime f2 CT < 1. 

The value of fl a is to be calculated at the decay epoch 
H ~ r. It is sometimes convenient to write 



T = 7 m 3 /A^P- 



(141) 



Then 7 ~ 10 -2 corresponds to gravitational-strength de- 
cay [94( and one expects 7 > 10~ 2 . 

Suppose first that there is continuous radiation domi- 
nation during the oscillation. If fl a <C 1, [27| 



Mpmj 1 / 2 

An approximation valid for any f2 CT is therefore 



(142) 



C 2 

M p to7 1/2 = Ml\l — 



Requiring a decay rate of at least gravitational 
strength, the first equality implies 



C 2 > lQ- l M P m. 



(145) 



Requiring that the decay takes place before the onset 
of nucleosynthesis, corresponding to p 1 / 4 > lMeV, the 
second equality implies 



C 2 > IQ-^M^/m 1 ' 2 . 



(146) 



11 According to the definitions made in this paper, £i n f = Cinfj but 



These bounds cross at m ~ 10 4 GeV, implying C > 
10 11 GeV. It will be important later that C might be 
either bigger or smaller than H. 

If pinf has a matter component C is bigger. In partic- 
ular, a contribution of curvaton particles, denoted by p c , 
gives 



C* 2 ~ Alp my 



1/2 _ 



M 2 n c , 



(147) 



where f2 c is evaluated at the onset of the oscillation. (The 
useful parameterization l|143fl was first given in p| , keep- 
ing just the contribution of curvaton particles.) 

Using these equations we arrive at the master formula; 



2fl a 



8 Or, 



u i 5o c 

q — 



(148) 



After adjusting the normalization of 5a this has the form 
Eq. (JTJ. Then the correlators are given by Eqs. (|14() - (|23|l . 



C. Special cases 

If we consider a single scale, Eqs. lf3"5Tl - (|4*T|) apply, and 
if in addition this scale is taken to be k ~ Hq(^ L^ 1 ) we 
can write things in terms of er; 



2Q a 



5a u f 5a 

q— + 7T «— 
a 2 \ a 



(149) 



with p given by Eq. I|135[) . With that understanding let 
us evaluate the correlators in some special cases. 

Let us assume that £ CT is dominated by the linear term, 
corresponding to 



,9V (H 

a 2 V 2tt 



< 1. 



( 143 ) Then 



(144) 



1/2 _ 2VLgq H 
^ 3p 2na 

3 3 P u /V Ca 

5 /NL ~ \~v7 



7"NL 



(36/25)/ 2 L 



(150) 

(151) 
(152) 
(153) 



In this case (, a can be the dominant contribution, de- 
manding e Nt " ~ 1. Then, unless u is extremely small, 
p ~ 1 leading to 



1/2 _ 2tt cr q H 

C ~~ 3 2lT(7 

3 3u 

t/nl — 



4fi ff 
(36/25)/ 2 L . 



(154) 

(155) 
(156) 
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For q = u = 1, corresponding to negligible evolution 
or evolution under a quadratic potential, Eqs. (|154f) and 
(I155|l reduce to the standard result |29| : 



,1/2 



Q.„H 



3na 
3, 3 1 

5 /NL _ ZnT' 



(157) 
(158) 



the bound on /nl requiring > 10 2 . 

Supposing further that the evolution actually is negli- 
gible, 



n a = 



a 2 + C 2 ' 



(159) 



This is the simplest version of the curvaton model. It 

gives H > T>l ,2 C, and then Eq. flU gives the 
bound 



D. 



H > 10 7 GeV. 



The case fi CT = 1 



(160) 



In the limiting case where f2 CT is indistinguishable 
from 1 there is no need of the sudden-decay approxima- 
tion. Before curvaton decay is appreciable, the curvaton- 
dominated cosmic fluid has P = making £, and hence 
C CT , a constant. The local value of a provides the ini- 
tial condition for the evolution of the separate universes, 
making them identical. As a result, £ remains constant 
throughout and after the curvaton decay process. 

To evaluate the non-gaussianity in this case one can 
use the SN formalism, which for the curvaton model is 
equivalent to using second-order cosmological perturba- 
tion theory. Adopting the small-tilt and p = 1 assump- 
tions, the calculation described in |3J applies so that 



r/l 



NL 



(161) 



With negligible evolution or a purely quadratic potential, 
/nl = — | ■ It would be interesting to know if such a small 
value will ever be observable. 



density on a slice of uniform energy density. Given the 
separate universe assumption, the inflaton perturbation 
5(j> cannot generate an isocurvature perturbation since 
it just corresponds to a shift back and forth along the 
inflaton trajectory. Any orthogonal light field Cj might 
create an isocurvature perturbation, and the same field 
might give the dominant contribution to the curvature 
perturbation so that the two perturbations would be fully 
correlated. This has been called a residual isocurvature 
perturbation |U EJ il] . 

In the curvaton model, a residual isocurvature pertur- 
bation obviously cannot be created after the curvaton 
decays. If the cdm or baryon number is is created by 
the curvaton decay an d p^ f contains no curvaton parti- 
cles, the argument of [23 gives a residual isocurvature 
perturbation S ~ —3(1 — fl<j)(- This is viable only if 
il a is close to 1. Repeating the argument of [2£j for the 
case that pi n { contains curvaton particles, one easily sees 
that they should be discounted when evaluating in 
the expression for S, allowing a true Sl a <C 1. 

Finally, suppose that the cdm or baryon number is 
created before curvaton decay. Then, whether or not p- m { 
contains curvaton particles, the argument of |2£j gives a 
residual isocurvature perturbation S/( ~ —3(1 — fi" ca ) 
where the superscript denotes the epoch of creation. 

The above formulas apply also to the fractional isocur- 
vature perturbation in the lepton number density, from 
which the neutrino isocurvature perturbation can be cal- 
culated |2^. It is not out of the question to generate 
the big lepton number density, that is needed to give a 
significant neutrino isocurvature perturbation [Sfij ]. 



F. The Bunch-Davies case 

So far a is unspecified. This is the unperturbed value 
of the curvaton field when the minimal box leaves the 
horizon. Now we consider the case that a 2 is equal to the 
variance (H/2tt) 2 jt a of the Bunch-Davies distribution. 
As we saw earlier, it becomes simpler in that case to use 
a maximal box, such that a = and 



= (^e> = 



2^ 



1 



(162) 



E. Induced isocurvature perturbations 

The status of isocurvature perturbations in the cur- 
vaton model is considered elsewhere (2^, |85| on the as- 
sumption that pi n f contains no curvaton particles. Let us 
briefly reconsider the situation when that assumption is 
relaxed. 

As defined by astronomers, an isocurvature perturba- 
tion S may be present in any or all of the baryon, cold 
dark matter or neutrino components of the cosmic fluid 
when cosmological scales first approach the horizon, be- 
ing the fractional perturbation in the relevant number 



Then the master formula can be written 

, (Sot) 



(163) 
(164) 



In the following I set equal to 1 the evolution factor qup. 
After adjusting the normalization of 5a, Eq. (|111J) gives 



V 1 / 2 = in a t\l 2 e~ m ° { ± 



k 

Ho 



(165) 
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If V{ = , observation requires Nt a < 1 and 
Eqs. CIT}-|TT3J) give 



pi/a 

3 f 
5 



2 »_ 



40„ 



Tnl = (36/25)/ NL 2 . 



(166) 

(167) 
(168) 



In this case t CT ~ 10 9 //nli which has to be very small 
indeed as was first noticed by Postma jf3 ■ 

If V C<T < Vt, strong tilt is allowed. Then Eqs. (fTT^) 
and l|113|) give 



5 /nl " ^ {~a 

tnl = 164 ( % 



-3iVt„ 



(169) 



-4JVt 



< 



The present bound /nl < 121 requires fio-e - ^*' 7 
9 x 1CU 6 and the present bound t^l < 10 4 requires the 
stronger bound ^l tT e~ Nta < 2 x 10~ 6 . These bounds are 
an extension of Eq. 1(54(1 , derived now for a maximal box. 
As in that case, they apply only on large cosmological 
scales, allowing a large curvature perturbation n much 
smaller scales. 

To summarize, we find in the Bunch-Davies case strong 
non-gaussianity on large scales (k ~ Hq) provided that 
Nt a > 1. This is because a typical region of size Hq 1 
then becomes very inhomogeneous by the end of relevant 
inflation. The non-gaussianity is reduced as the scale is 
decreased, and (with fl a = 1) is small on the scale leav- 
ing the horizon at the end of relevant inflation precisely 
because a typical region of size fc" 1 is still quite homoge- 
neous. 

Although the maximal box provides the neatest result, 
it is interesting also to see what happens with a minimal 
box. For a box of any size, Eq. I) 163(1 becomes 



1 



(171) 



Setting u = q = p = 1 and adjusting for the normaliza- 
tion of 5a, we can calculate the correlators of £ a from 
Eqs. 1(371) . 1(38(1 . arid (|40|l and then take their expectation 
values within a maximal box, to find 





= 4t CT ( 






= K ( 




(tnl) 


= 164 





-3JVi„ 



/ (172) 



A] •r-'/CI7:-!) 



2/„ 



-4JVt„ 



/ 



fc 

(I/fn) / ( -jj- 



3/„ 



^- 3 il74) 



(175) 



In accordance with the discussion of Section III Dl these 
expressions are independent of the box size. But the 



split into the linear plus quadratic term depends on the 
box size. With a maximal box the linear term vanishes, 
even in gaussian regime Nt a < 1. With a minimal box 
the linear term dominates on large scales, even in the 
non-gaussian regime. 

With a minimal box, negligible tilt and C = Co-) 
Eq. (|157|l applies. Inserting the Bunch-Davies expecta- 
tion value for a 2 then reproduces Eqs. (|166fl - l|168|) . If in- 
stead we consider the PNGB case, assuming A <C H* so 
that a has a flat distribution up to a max = /, Eq. 1)166(1 
becomes 



^/ 2 >^4i4 >>f7CT A 



A / 

Using Eq. i(119|) we have again ^l rT ^ 
indistinguishable from 1. 



/ 



(176) 



< 1, making t a 



G. Cosmic uncertainty 

To make contact with previous work, I assume in this 
subsection that the curvaton field has negligible spectral 
tilt, and has negligible evolution before the oscillation 

1 /2 

starts. Also I consider Q = (2/5)7-^ whose observed 
value is 2 x 10~ 5 . 

We are not going to be concerned with precise values, 
but in this context we do not want to exclude the case 
a 2 < jj2 rp Q handle it W e can replace a in Eq. 1(157(1 by 
\J a 2 + H 2 , leading to 



Q a 2 + H 2 + C 2H n ° 
a 2 + H 2 



H 



y/a 2 + H 2 



NL 



a 



a 2 + H 2 + C 2 



(177) 



(178) 



In this expression, H is evaluated during inflation. 

In the Bunch-Davies case, the probability distribution 
for a is gaussian, making it more or less flat up to a 
maximum value of order the variance; cr max ~ H/ \/t^ 3> 
H. In the PNGB case (Eq. (fTT£|l with A < H) the 
the probability distribution is almost perfectly flat, up 
to a smaller maximum which could be much less than H. 
Using the flat distribution, one can work out the non-flat 
distribution for the correlators. To discuss this, I take 
the tilt t ff to have a small fixed value, consistent with 
observation. 

There are three parameters C, H and a and I will 
divide the parameter space into regions, separated by 
strong inequalities to allow simple estimates. 

Consider first the case C <C H. In the regime H <C 
a <C cr max we have Q ~ H/a, with fl a very close to 1 
and /nl = —5/4. In the regime a -C H, the quadratic 
term dominates C making Q ~ 1 in contradiction with 
observation. Given the flat distribution for a, the lengths 
Cmax — H and H over which a runs for the two cases give 
their relative probabilities for a randomly-located region. 

Now consider the case C H. There are three 
regimes; 
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1. The nearly gaussian regime C <C o < <7 max . Here 
Q ~ if/cr and fl a ~ 1 making / NL — —5/4. 

2. The strongly non-gaussian regime < c C C. 
Here Q ~ iJa/C 2 and ~ CT ~ ct 2 /C 2 < 
1. The lower part of this range is forbidden by 
observation. 

3. The regime ct <C H. Here the quadratic term of 
Eq. I|148|) dominates leading to Q ~ if 2 /C 2 , and 
to /nl ~ 1/Q in contradiction with observation. 

The relative probabilities, that a given region corre- 
sponds to one or other of these cases, are proportional 
to the intervals Act given at the beginning of each item. 
Within each case, the probability distribution for Q is 

dP = da = ^rdQ. (179) 
dQ 

H. Anthropic considerations 

So far we did not take into account anthropic consider- 
ations. They have gained force recently, since it appears 
that string theory allows a very large number of field the- 
ory lagrangians. Both the idea and the methodology of 
anthropic arguments are very controversial, as empha- 
sized for instance in §, but let us proceed. 

Anthropic arguments suggest [95| that Q has to be 
in a range Q m i n < Q < Q max corresponding roughly 
to 1CP 6 < Q < 10~ 4 . If the cosmological constant is 
taken to be fixed the probability distribution within this 
range is more or less flat. But Weinberg argued [9(| that 
the cosmological constant itself should be regarded as 
having a flat probability distribution, since there appears 
to be no theoretical argument that would give a definite 
value, in particular zero. He showed, before the data 
demanded it, that anthropic arguments suggest a value 
appreciably different from zero. As summarized in 
subsequent studies have shown the preferred value to be 
compatible with observation. Accepting this viewpoint 
for the cosmological constant, it has been argued [o3| 
that the probability distribution of Q within the above 
range is 

dP oc Q 3 dP prim , (180) 

where the prior dP pr i OI is the probability distribution if 
anthropic considerations are ignored. I will use this esti- 
mate in the following discussion, without tryin g to take 
on board the impact of some more recent work |98( . 

At this point, one may wonder why the focus is exclu- 
sively on the overall normalization Q 2 of the spectrum. 
What about the spectral tilt, and the measures /nl and 
tnl of non-gaussianity? The original arguments of Harri- 
son and Zeldovich 100] for tilt in the range \t a \ < 1 
may be regarded as anthropic, but we are now dealing 
with an observational bound more like \t a \ < 0.01. It 
seems clear that no anthropic consideration will directly 



produce this result, and there was no objection to val- 
ues \t a \ ~ 1 before the cmb anisotropy ruled them out. 
Coming to non-gaussianity, it is again hard to see how an- 
thropic consideration will directly constraint it, and there 

— 1/2 

was no objection even to the extreme case /nl ~ 
until it was ruled out by observation. From the anthropic 
viewpoint, the small tilt and non-gaussianity presumably 
are produced accidentally, by anthropic constraints on 
other parameters including Q. 

If Q is generated by the inflaton perturbation in a 
single-component model, it depends almost entirely on 
some parameters in the field-theory lagrangian. (There 
is some dependence on the post-inflationary history via 
N but it would take a big variation of that history to 
have much effect on N.) Then the tilt also depends only 
on the field theory parameters, while the non-gaussianity 
is automatically negligible. If the field theory parameters 
were taken to be fixed, the prior would be a delta function 
and there would be no room for anthropic arguments. 

In contrast, curvaton-type models depend also on the 
background values of one or more fields, and if inflation 
begins early enough one has no option but to consider 
their variation within the very large and smooth inflated 
patch that we occupy. This was pointed out some time 
ago in |63,|73 for the PNBG case, and has been discussed 
more recently in [j| for the Bunch-Davies case. The point 
here is that anthropic considerations concerning Q may 
demand, or anyhow favour, a value ct far below cr max . 

A precisely similar situation exists with respect to the 
nature of the CDM. If it consists of neutralinos, or of ax- 
ions created by the oscillation of cosmic strings, the CDM 
density is given in terms of parameters of the lagrangian. 
If instead it is the oscillation of a nearly homogeneous ax- 
ion field which existed during inflation, the CDM density 
varies with our location within the inflated patch (and 
so does the CDM isocurvature pe rturb ation which is in- 
evitable in that case 0). Linde jlOlj provided the first 
concrete realization of anthropic ideas, when he pointed 
out that we might need to live in a place where the ax- 
ion density is untypically small. The anthropic probabil- 
ity for the axion density has recently been investigated 

Now I analyse the situation for the actual curvaton 
model, generalizing two recent discussions 0, 0. The 
spectrum Q 2 is given by Eq. (|177f) . and we are assuming 
that the probability distribution is flat within a range 
< ct < CT max . There is also the the anthropic constraint 
Qmin < Q < Qmax- These inequalities define a rectangle 
in the Q-a plane, and we can only use the part of the 
curve l|177|) that lies within this rectangle. 

The location of the rectangle relative to the curve de- 



In this case the axion density is proportional to (o -2 + H 2 ) where 
<T is the average axion field in our part of t he U niverse and H 2 
com es from the long wavelength fluctuations JrjJ . The authors of 
98] drop the H 2 term, which might possibly affect their results. 
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pends on the parameters C and H which define the curve, 
and on cr max whose value was discussed earlier. In this 3- 
parameter space, there will be an unviable regime where 
no part of the curve lies within the rectangle. In the 
opposite case, part of the curve is within the allowed 
rectangle, and putting Eq. Ij 179(1 into Eq. (|180J) gives the 
probability distribution for Q; 



dP cx Q 3 (a)da cx Q 3 ^-dQ 



(181) 



Consider first the case C <C H. The regime a < H 
is unviable because it gives Q ~ 1 > Q max - Therefore 
we consider the regime a 3> H, and assume that cr max is 
big enough that its value does not matter. (As <7 max is 
reduced from some such value, it remains irrelevant to the 
following discussion until the allowed rectangle ceases to 
intersect the curve making the model unviable.) This is 
the version of the curvaton model considered by Garriga 
and Vilenkin Q. As we noted already it corresponds to 
Q ~ H/a with £l a very close to 1 and /nl = —5/4. As 
a is reduced, Q rises to a maximum of order 1, but only 
the regime Q < Q max is allowed. Within this regime the 
probability is 



dP (x QdQ, 



(182) 



making Q — Q max the most likely value. Taking Q max = 
10~ 4 , the probability that Q is at or below its observed 
value is 1/25. As these authors emphasize, the low prob- 
ability for the observed value need not be taken very 
seriously if only because Q m ax is not at all well-defined. 

Let us move on to the regime C ^ H . In this case 
Q(a) has a maximum at a ~ C, with the value H/C. 
Let us consider first the case that cr max is big enough for 
its value not to matter. 

If H/C < Q max , the peak value is anthropically al- 
lowed, and because of the Q 3 factor is in fact favoured; it 
had better agree with observation if the anthropic argu- 
ment is to work. It corresponds to f2 ff ~ 1 and /nl ~ 1- 
This is the case considered by Linde and Mukhanov 
Now suppose instead that H/C is bigger than Qmax, so 
that a region around the peak is excluded. The proba- 
bilities to the right and to the left of the peak are 



(iPright 

dPich 



H 3 



-da 

o 

Ha 
~~CP- 



da. 



(183) 
(184) 



Integrating these expressions, the the relative probability 
for being in the two regions are 



flcft 
-fright 



H 



< 1. 



(185) 



The right-hand region is therefore preferred anthropi- 
cally, leading again to the estimate l|182fl . 

We have still to consider the case that a < a ma x is 
a significant constraint. As <7 m ax moves down from a 



large value it will at some point exclude the right hand 
part of the curve. Then Eq. I|184|) applies which gives 
dP cx Q 3 dQ. The probability that Q is at or below the 
observed value is now only 1.6 x 10~ 3 (with Q max = 1CP 4 ) 
which might be regarded as a catastrophe for anthropic 
considerations. And we are now in the regime of strong 
non-gaussianity, which means that depending on param- 
eters there might be a violation of the present observa- 
tional bound on /nl- 

As <r max moves further down, it will start to cut into 
the left hand part of the curve. Eventually, the peaked 
probability distribution that we encountered in the pre- 
vious paragraph is CUt off at Qmax = Q(c m ax) < Qmax so 

that Qmax becomes the preferred value, which had bet- 
ter agree with observation if the anthropic argument is to 
work. Again, one has to check that /nl is small enough. 

This completes our discussion of the anthropic status 
of a simple version of curvaton model. We see that the 
situation is rather complicated. It would get still more 
complicated if we allowed evolution of the curvaton (not 
to mention the possibility that the curvaton contribution 
is sub-dominant) or if we considered a very non-flat prior 
probability for a such as might come a departure from 
the probability distribution (|117fl . 



I. Comparison with a previous work 

To a considerable extent Sections llV EH1V fll represent 
a development of p| (see also |2(|). Where they are com- 
parable, the results are in broad agreement. 

An expression essentially equivalent to Eq. i|165|) is de- 
rived in [2(j. To be precise, the expressions formally 
coincide because our factor e Nt " is equal to their fac- 
tor (H^/Xo) defined in @. But we make a distinction 
between the number of e-folds of inflation (after the ob- 
servable Universe leaves the horizon) and the number of 
relevant e-folds, because it is the latter that should be 
identified with N. 

The only other significant difference is one of interpre- 
tation, concerning the Bunch-Davies gaussian field cr(x) 
which they call the curvaton web. One's view about the 
curvaton web depends on the interpretation of a. Within 
the observable Universe, the local value of a will vary 
from place to place, and its variation may be observable. 
The variation of a might be important if, for instance, 
one is considering the galaxy distribution in a relatively 
small region surrounding our galaxy. Indeed, the av- 
erage of a within this region (at horizon exit) may be 
significantly different from its average within the whole 
observable Universe. (Analogous considerations for the 
axion dark matter density were pointed out for instance 
in 17].) With this interpretation of a, the steep spatial 
gradient of <r(x) evident in one of the simulations may 
be an observable effect, as the authors remark. 

On the other hand, the analysis of local effects within 
the observable Universe is a rather tricky business even 
if one is not concerned with a varying scalar field (but in- 
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stead with say the local expansion rate) . It may therefore 
be useful to have a division of labour, whereby models 
of the early universe give the correlators defined with 
the box size comfortably enclosing the whole observable 
Universe. These then provide the starting point for the 
analysis of local effects, which can be done at a later and 
more sophisticated stage of the research. Certainly that 
is the viewpoint usually taken when the curvature per- 
turbation is supposed to come from the inflaton, and it 
has been the viewpoint also of the present paper. From 
this viewpoint, observation is sensitive to just one point 
on the curvaton web, whose spatial gradient ceases to be 
physically significant. The gaussian probability distri- 
bution for a within the minimal box, when inserted into 
Eqs. ifTTTjl and lfT?%|) directly gives the cosmic uncertainty 
of the correlators, and with the usual 'cosmic variance' of 
the almost-gaussian CMB multipoles this covers all pos- 
sibilities for what will be observed even though it may 
take some effort to work them out. 



V. CURVATON- TYPE MODELS AFTER WMAP 
YEAR THREE 

If a curvaton-type contribution Q a dominates the cur- 
vature perturbation, the t ensor fraction is tiny. Then the 
WMAP year three results ^02j combined with the SDSS 
galaxy survey give n — 1 ~ — 0.0521q'q]J, and the result 
hardly changes if WMAP data are used alone or with 
several other relevant data sets. 

If this measurement of a small negative spectral tilt 
holds up it has important consequences for curvaton-type 
models for the origin of the curvature perturbatio n. A s 
was noticed in the early days of their exploration |103| , 
the most natural expectation for these models is that the 
spectral tilt t a of the curvaton-type field a is practically 
zero. This is because, in contrast with the inflaton poten- 
tial in a non-hybrid model of inflation, the potential of a 
does not 'know' about the end of inflation. The potential 
already has to be exceptionally flat just to convert the 
vacuum fluctuation of a into a classical perturbation, and 
in the absence of any reason to the contrary one might 
expect that the departure from flatness will be too small 
to observe. 

If this expectation is accepted, the spectral tilt pre- 
dicted by a curvaton-type model is 

n — 1 = — 2e* curvtiltO , (186) 

where e is the flatness parameter of slow-roll inflation, 
or more generally is the parameter e# = —H/H 2 . To 
reproduce the observed tilt we need a more or less scale- 
independent value e ~ 0.025 

A mong a suite of models considered in a recent sur- 
vey |104| . the large-field models with V oc cf> a (chaotic 
inflation) give roughly the correct e. The degree of tilt 
depends on N, defined in this context as the number of 
e-folds of inflation after the observable Universe leaves 
the horizon. The best-motivated case is a = 2, because 



it may be obtained as an approximation to Natural In- 
flation. Taking N = 50, this gives n — 1 = —0.020 which 
is a bit too small compared with the observed value. The 
multi-component version of this potential can help by re- 
ducing Ti — l, but no investigation has been done to see 
how far one can go in that direction. 

Increasing a increases e by a factor a/ 2, so that a — 4 
or 6 give a tilt agreeing with observation. (If the inflaton 
perturbation is required to generate the curvature per- 
turbation such values of a are excluded, but that is not 
the case here.) The problem with increasing a is that 
it lacks motivation, either from stri ng th eory or from re- 
ceived wisdom about field theory |l04j . If one accepts 
an increased a it may be more sensible to regard a as a 
non-integer, providing an approximation to the potential 
over the relevant range of 4>. 

What about the possibility of giving the curvaton-like 
field a a significant negative tilt t ai say enough to allow 
?i — 1 ~ t a ? This requires a to be on a concave-downward 
part of its potential during inflation. Taking the view 
that the value of a is an initial condition to be assigned 
at will (possibly with anthropic restrictions) this need 
not be a problem. If instead the value has the stochastic 
probability distribution l|117(l , there is some tension as we 
saw after Eq. (|118fl but one can still achieve the required 
small tilt with reasonable probability. Only in the sim- 
plest version of the actual curvaton model does Eq. I|117|) 
(with a quadratic or periodic potential) demand negligi- 
ble tilt t a . In general then, curvaton-type models can 
easily give the curvaton perturbation a suitable negative 
tilt allowing n — 1 ~ t a . The challenge in that case 
though, is to explain the actual value n — 1 ~ —0.05 in 
a natural way. As was noticed a long time ago [||, a 
wide class of inflation models generating the curvature 
perturbation from the inflaton do just that, by making 

77-1 1/N. 

If curvaton-type models are rejected as the origin of the 
curvature perturbation, that rejection is itself a powerful 
constraint on any early-universe scenario involving scalar 
fields other than the inflaton. Such fields must either be 
heavy during inflation, or else their contribution to the 
curvature perturbation must be negligible. The possible 
problem caused by a curvaton-type contribution being 
too big is analogous to the moduli problem, and indeed 
might even be part of that problem [2^, . 



VI. CONCLUSION 

The main result of Section [H] is the extension of Q 
to include strong spectral tilt, which is allowed if the 
curvaton contribution gives a sub-dominant contribution 
to the curvature perturbation. 

Section II I II shows how to derive Eq. JTJ and its gen- 
eralizations, in the presence of both spectral tilt and 
non-gaussianity. We note that the curvature perturba- 
tion generated by curvaton-type model long after infla- 
tion might have a very low ultra-violet cutoff fc max . This 
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means that the scale of inhomogeneities in a matter com- 
ponent of the cosmic fluid might not extend much be- 
low the scale 10 6 M Q required for the formation of the 
first baryonic objects. On the other hand it might, in 
which case the curvaton-type contribution to £ might be 
negligible on large cosmological scales, but big on sub- 
cosmological scales even allowing primordial black hole 
formation. 

On the technical side, we note that the horizon- 
crossing formalism does not work when spectral tilt and 
non-gaussianity are both to be included. Instead one 
should evaluate the field perturbations at a fixed epoch, 
which might as well be the one when /e max leaves the 
horizon. 

Section IIVI revisits the actual curvaton model, taking 
into account all of the possible effects that have been 
noticed. If the curvaton contribution to the curvature 
perturbation has strong positive tilt, it can be negligible 
on cosmological scales but big enough to form primordial 



black holes on smaller scales. Recent discussions concern- 
ing cosmic uncertainty and the anthropic status of the 
curvaton models are extended. This section also demon- 
strates that the prediction /nl = —5/4, of the simplest 
version of the curvaton model, can be obtained with- 
out recourse to the sudden-decay approximation. Conse- 
quently, a detection /nl = — 5/4 would be a smoking gun 
for the simplest version of the curvaton model. Finally, 
in Section we looked at the status of curvaton-type, in 
the light of the recent measurement of negative spectral 
tilt for the curvature perturbation. 
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